Zeeman and Orbital Effects of an in-Plane Magnetic Field in Cuprate Superconductors 

Kun Yang 
National High Magnetic Field Laboratory and Department of Physics, Florida State University, Tallahassee, Florida 32310 






S. L. Sondhi 
Department of Physics, Princeton University, Princeton, New Jersey 08544 

(February 1, 2008) 



o 

o 
o 

u 

Oh' 



X3 

o 



> 

OO 

o 

0^ 
OS 
-I— > 



I 

O 

o 



X 



We discuss the effects of a magnetic field applied parallel to the Cu-0 (ab) plane of the high Tc 
cuprate superconductors. After briefly reviewing the Zeeman effect of the field, we study the orbital 
effects, using the Lawrence-Doniach model for layered superconductors as a guide to the physics. 
We argue that the orbital effect is qualitatively different for in-plane and inter-layer mechanisms for 
superconductivity. In the case of in-plane mechanisms, interlayer couplings may be modeled as a 
weak interlayer Josephson coupling, whose effects disappear as H —> oo; in this case Zeeman domi- 
nates the effect of the field. In contrast, in the inter-layer mechanism the Josephson coupling is the 
driving force of superconductivity, and we argue that the in-plane field suppresses superconductivity 
and provides an upper bound for Hc2 which we estimate very crudely. 



One of the most important milestones in the study of 
the cuprate superconductors is the identification of the 
predominantly 6.^2 _y2 symmetry of their pairing order 
parameters. This discovery has motivated a tremendous 
amount of theoretical work on the physical properties of 
superconductors with unconventional pairing symmetry. 

Among them, it was recently pointed outEJ that the re- 
sponse of a two-dimensional (2D) dx2-y2 superconductor 
to the Zeeman coupling between the spins of the elec- 
trons and an external magnetic field is much stronger 
than that in an s-wave superconductor. This is so, for 
unlike the fully gapped s-wave case, in a dx2-y2 su- 
perconductor there exists nodal points in the order pa- 
rameter and therefore gapless quasiparticle excitations. 
Consequently, an arbitrarily weak Zeeman field gener- 
ates a finite density of spin-polarized quasiparticles in the 
ground stateu. It was also shown that the Zeeman field 
enhances the low-T specific heat, thermal conductivity, 
and tunneling deasity of states, while suppressing the 
superfluid densityld'El. At strong fields, the system sup- 
ports a FuldCpRerrell-Larkin-Ovchinnikov type supercon- 
ducting state.l3'EJ Eventually, a sufficiently strong Zeeman 
field suppresses superconductivity completely, leading to 
a "Pauli limit" on the upper critical field. 

The results mentioned above are expected to be partic- 
ularly relevant to the cuprate superconductors when the 
external field is oriented parallel to the Cu-0 (ab) planes 
(for possible caveats on the influence of disorder see [ ^), 
as the cuprates are quasi-2D systems with electrons pre- 
dominantly moving in the Cu-0 planes. In a truly 2D 
system, the orbital motions of the electrons would not 
sense the existence of a parallel magnetic field, and the 
Zeeman effect would be the only effect of the field. Ex- 
perimentally, there is some evidenceO suggesting that the 
upper critical field for the parallel orientation is indeed 
Pauli limited in YBCO at low temperatures. 

However for smaller fields there do exist orbital effects. 



even when the field is parallel to the planes, due to the 
presence of inter-plane couplings. It is the purpose of 
this paper to report results of studies of such orbital ef- 
fects, and to compare them with the Zeeman effect in 
field ranges that are relevant to current experiments. We 
will argue that the nature and importance of the orbital 
effects depend sensitively on whether superconductivity 
is driven by in-plane correlations or inter-plane hopping. 
Thus their study may provide a way to distinguish be- 
tween in-plane and inter-plane mechanisms for supercon- 
ductivity in the cuprates. 

The orbital effects of the parallel field may be sum- 
marized qualitatively as the following. The field gener- 
ates screening currents both within and across the planes. 
The current within the plane|-generates a "Doppler shift" 
of the quasiparticle spectrumEl, and therefore a finite den- 
sity of states (DOS) for the quasiparticles. This current- 
induced DOS is proportional to VH for -ff ||a3, and thus 
dominates the Zeeman effect (which gives rise to a quasi- 
particle DOS that is proportional to Ba) at low H. We 
show below that the situation is quite different for an in- 
plane field (iJ_Lc). Also the inter-plane Josephson cur- 
rent suppresses the inter-plane Josephson coupling en- 
ergy. Thus, depending on how crucial the Josephson cou- 
pling energy is to superconductivity itself, the response 
of the system to the in-plane field can be very different. 

We [-.model the cuprates by the Lawrence-Doniach 
modeB that is appropriate for layered superconductortfl: 
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Here n is the layer index, a oc T — Tc, and other param- 
eters may be assumed to be T independent. For T very 
close to Tc, the coherence length becomes very large, so 



that continuum approximation can be used along the c 
direction, and (|l|) reduces to the usual Ginsburg-Landau 
theory with anisotropic mass, and the mass along the c 
direction is rric = ^^ where s is the spacing between 
layers. The parameter 7^ = rric/mab is a measure of the 
degree of anisotropy of the system; for examples, 7^ ~ 50 
for YBCO, and 7^ « 20000 for BSCCOEI. 

We are primarily interested in the low T regime, where 
the layered structure is important; more precisely we are 
interested in the case T < T* , the latter being the tem- 
perature at which the orbital Hc2 diverges. At such 
temperatures we distinguish two regimes: a low field 
regime in which the current flowing in the layers in- 
creases with H on account of increasing vortex density, 
and a high field regime in which the increasing density 
is counteracted by the cancelation between neighboring 
vortex uatterns. The former has also been discussed by 
VolovikD and also yields a DOS ex \/H. Here we compute 
the answer at high fields. We assume the external field 
H = H{1,Q,0) is along the a direction. It is known that 
Hci is very low (< lOOGauss) for a parallel field, and 
for fields of order H ^ IT (which is the range of interest 
in the present paper and typical for experimental stud- 
ies), the field basically penetrates the system uniformly 
and the screening effect is negligible. (This justifies ne- 
glecting the field energy term in (|l|)). Thus we can as- 
sume the field inside the superconductor is the same as 
that of the external field and uniform. We use the gauge 
A„ = ns7J(0,I,0). 

We need to solve the equation for -^n: 
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For in-plane mechanisms, the Josephson coupling t is 
weak and secondary to superconductivity. To find the 
solution to (ni, we first set i = 0, and then treat t as a 
weak perturbation. For i = 0, all the layers decouple, 
and we obtain 



V'„(x,y) = v/^^e'(*"+'="^), 
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and (j)n are arbitrary phases. Upon 



introducing t these solutions spread to more than one 
layer aruilock phases, leading to the unique, leading order 
solutionllil 
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where Ak — 2esH/hc, £ ~ ^Jfic/2eH is the magnetic 
length for Cooper pairs, ?]„ — 1 for n — 4j or 4j' + 1, and 
— 1 otherwise (j is an integer). In order for the pertur- 
bative approach to be valid, we need to have 
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which is always valid in the strong field limit. 
The typical in-plane superfluid velocity is thus 
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and the Doppler shift energy isEHD 
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Since the Zeeman energy Ez — fJ-sH oc _ff, it always 
dominates the Doppler shift at sufficiently strong field. 

At field H = lOT, using ep ^ 2eV, kp ^ lA"^ and 
s ~ lOA, we obtain for BSCCO (in which (0) is satis- 
fied) Ed/Ez ^ 0.4, indicating Zeeman effect starts to 
dominate the orbital effect in this field range in BSCCO. 
(g) is not satisfied in YBCO at this field due to less 
anisotropy; but using the above results anyway (which 
results in an overestimate of the in-plane current) we 
find Ed/Ez ~ 150, suggesting that the orbital effect 
will dominate. 

In the above discussion we have assumed that the in- 
terlayer Josephson couplings are secondary to supercon- 
ductivity in cuprates, and may be treated as weak per- 
turbations. This is appropriate for in-plane mechanisms 
for superconductivity in cuprates—Hpwever it was pro- 
posed by Anderson and coworkerstSEJ that the interlayer 
Josephson coupling (or pair-hopping) may actually be 
the driving force behind superconductivity in cuprates. 
In this case, the interlayer Josephson coupling can no 
longer be treated as a weak perturbation as above, and 
we expect an in-plane magnetic field to have a much more 
dramatic orbital effect. This is because the in-plane field 
tends to suppress interlayer phase coherence, and there- 
fore Josephson energy and superconductivity in this case. 
Among other things, the orbital effect of the field pro- 
vides a mechanism for the upper critical field Hc2 ■ While 
for in-plane mechanisms, the field simply decouples 
layers, and ifc2 can only come from the ZejEaaan effectlH 

In the inter-layer pair hopping modeO'Ej, the multi- 
layer structure of the cuprates is important. In particu- 
lar, the Josephson coupling among layers within the same 
unit cell might be much stronger than between different 
cells, thus the large anisotropy could largely reflect the 
weakness of the inter-unit cell Josephson coupling only. 
For simplicity, we study in the following in a single bi- 
layer (as in Ref. |l^), and neglect the coupling between 
this pair of layers and the rest of the system. A proper 
treatment of even this problem requires a microscopic 
theory that does not really exist. So here we will content 
ourselves with a speculative illustration of what such a 
treatment might produce. To this end we again use the 
Lawrence-Doniach model, in a slightly different form and 
with a crucial reinterpretation of the parameters: 
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layer coherence, as the Josephson energy is the origin of 
superconductivity. Thus the solution in this case is 
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Here a' — a + 1. We have used a different gauge: For a 
parallel field H = Hy, we choose A = Hxz, so the A af- 
fects the phase of the interlayer Josephson coupling, but 
not in-plane kinetic energy. The physics, of course, is un- 
affected by the choice of gauge. The reason to isolate the 
last term out is that this term describes interlayer cou- 
pling only; while all the in-plane properties are described 
by other terms. 

For in-plane mechanisms, each individual layer is su- 
perconducting (at low T) without interlayer coupling. 
Thus we must have a' < 0, which means the reason that 
a superconducting order parameter ip is developed is be- 
cause of the lowering of in-plane free energy. Without 
a parallel field, the free energy can be further lowered 
by having the same phase for tpi and ^^2 (and uniform 
throughout the system so that there is no in-plane ki- 
netic energy cost), so that the systems gains Josephson 
coupling energy as well. 

In pair hopping mechanism, however, the individual 
layers are not superconducting without interlayer Joseph- 
son coupling. This means a' > 0! Instead the system 
becomes superconducting at if = by developing a non- 
zero ip in both layers, and making them phase coherent 
("01 — ip2), thereby gaining Josephson coupling energy. 
Thus as long as t > a', the system is in the supercon- 
ducting phase, despite a' > 0. 

As long as H = 0, the ground state of (pt) is the same: 
"01 = "02 = -\/(i — q')//3, and there is no obvious distinc- 
tion between these two cases; the system is phase coher- 
ent both within each individual layer and between the 
two layers. The situation is very different when H ^ 
and is sufficiently large. This is because for a loop that 
encloses a finite amount of flux, one must accumulate a 
finite amount of gauge-invariant phase difference propor- 
tional to the flux (see Fig. la). It is therefore no-longer 
possible to maintairuip-plane and interlayer phase coher- 
ence simultaneouslytJ. 

For a' < (in-plane mechanism), the situation is 
basically tlie same as a Josephson junction in a mag- 
netic fieldt2l; we have V'j = e*"^^ ^/—a'/P, with no cor- 
relation between (pi and (p2- The leading order cor- 
rection to this is proportional to t/H. The field in- 
duces an oscillatory Josephson current between the lay- 
ers: jz « {—a't/0)sm{^^j^+<pi~<p2), while the in-plane 
current goes to zero as 1/H at large H (see Fig. lb). 
The Josephson coupling term in (ra) averages to zero in 
this solution. In short, the system gives up interlayer co- 
herence (as manifested by the loss of Josephson coupling 
energy and appearance of Josephson current), in order to 
maintain in-plane coherence (therefore no in-plane cur- 
rent and kinetic energy). 

For a' > (interlayer mechanism), on the other hand, 
the system will do everything possible to maintain inter- 
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In this solution the interlayer coherence is maintained 
and there in no interlayer Josephson current induced by 
the field. As a price, the system loses in-plane coherence, 
and there is diamagnetic current (ex H) and in-plane ki- 
netic energy loss (ex H^) (see Fig. Ic). Clearly, as H 
increases, at some point the loss of in-plane kinetic en- 
ergy will overwhelm the gain of Josephson energy, and 
the system will cease to superconduct. This is the or- 
bital upper critical field H°2 (to be distinguished from 
the Zeeman upper critical field H^2)i ^s determined by 
the point where ip vanishes in (|9[). For in-plane mech- 
anisms, H°2 —> c». If we assume that interlayer mech- 
anism is responsible for superconductivity in cuprates, 
a crude estimate yields in BSCC02212 H°2 ^?pOT at 
T = 0. Unfortunately this is comparable to Tf^tlH. How- 
ever, in principle, the Zeeman effect can be suppressed 
by introducing spin-orbital scatterers to the system which 
could enable the two to be distinguished. 

The qualitative difference in the orbital response to a 
parallel magnetic field can provide a way to experimen- 
tally distinguish between in-plane and interlayer mech- 
anisms for cuprate superconductivity. We note the in- 
terpretation of previous experimental attempts to distin- 
guish between tiie two based on c-axis penetration depth 
is controversially. 
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FIG. 1. (a) A closed path that encloses a finite amount 
of flux, (b) The current pattern in a bilayer system in the 
presence of strong parallel magnetic field, for in-plane mech- 
anism of superconductivity, (c) Same as (b), for inter- layer 
pair hopping mechanism of superconductivity. 



